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Superconducting energy gap distribution in c-axis oriented MgB2 thin 

film from point contact study 
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We have analyzed about hundred voltage-dependent differential resistance dV/dI(V) curves of metallic 
point contacts between c-axis oriented MgB2 thin film and Ag, which exhibit clear Andreev reflection fea- 
tures connected with the superconducting gap. About one half of the curves show the presence of a second 
larger gap. The histogram of the double gap distribution reveals distinct maxima at 2.4 and 7 meV, while 
curves with a single-gap features result in more broad maximum at 3.5 meV. The double-gap distribution is 
in qualitative agreement with the distribution of gap values over the Fermi surface calculated by H. J. Choi 
et al. (cond-mat/0111183). The data unequivocally show the presence of two gaps As = 2.45 ± 0.15 meV 
and Al = 7.0 ± 0.45 meV in MgP>2 with gap ratio Al/As ~ 2.85 ± 0.15. Our observations prove further 
a widely discussed multigap scenario for MgB2, when two distinct gaps are seen in the clean limit, where a 
single averaged gap is present in the dirty one. 

PACS: 73.40.Jn, 74.76Db, 74.80.Fp 



Introduction. Direct spectroscopic investigations of 
the superconducting order parameter in recently dis- 
covered [0 superconductor MgB 2 with T c ~40K by 
tunneling P, |s|, f|, p, 3|, [7, B, p|, and point-contact 
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spectroscopy [jlO, [H, l|, p| 
ambiguously an energy gap A in the quasiparticle den- 
sity of states (DOS). However, the experimental re- 
sults are controversial as to the gap width A, whose 
variation from 1.5 to 8 meV (see, e. g., review Jl8|) 
is unexpectedly large, pointing to the possibility of 
multiphase or nonhomogeneous samples, degraded sur- 
face, or anisotropic energy gap. Another way to solve 
this puzzle is to consider two superconducting gaps 
■in MgB 2 , as proposed by Liu et al. ]19| , accounting 
complex electronic structure of MgB2 with both quasi- 
2D and 3D Fermi surface f20|| . Indeed, several papers 
§, @, §, have reported double gap struc- 

ture in the differential conductance (resistance) with the 
smaller gap being far below weak-coupling BCS value 
A=1.76keT c ~6meV and the larger gap slightly above 
the standard BCS one, in accordance with theory [fL9[ - 
Therefore, one of the intriguing key issues of super- 
conducting state of MgB2 is whether the double gap 
state is intrinsic or the spread of the gap values is a re- 
sult of anisotropy, nonhomogeneity, surface effect, etc. 
In other words, before macroscopic high quality single 
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crystals will be available for thorough investigations, the 
sample imperfection may raise doubts about the final 
conclusion. However, in our mind, good reproducibil- 
ity of the double-gap values given by different authors 
§, §, ||, HI, by different, in their physical back- 
ground, methods such as tunneling and point-contact 
spectroscopy carried out on different samples such as 
pellets, films, grains, all this with a great degree of prob- 
ability supports intrinsic nature of the double gap in 
MgB 2 . 

In this paper we will give further confirmation of 
double gap scenario in MgB 2 based on analysis of about 
hundred point-contact spectra of c-axis oriented thin 
films. 

Experimental and calculation details. We have mea- 
sured the high-quality c-axis oriented 0.4 /im thick 
MgB 2 film || grown by a PLD technique on Al 2 03 
substrate. The resistivity of the film exhibits a sharp 
transition at 39 K with a width of ~ 0.2 K from 90% 
to 10% of the normal state resistivity pl[ . The residual 
resistivity p Q at 40 K is - 6 ^cm0 and RRR=2.3. 

Different point contacts (PCs) were established in 
situ directly in liquid 4 He by touching as-prepared sur- 
face (sometimes etched by 1% HC1 solution in ethanol) 
of the MgB 2 film by a sharpened edge of an Ag counter- 
electrode, which were cleaned by chemical polishing in 

2 ) There is a scattering by factor of 4 in po f° r the different films. 
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HNO3. This geometry corresponds to the current flow- 
ing preferably along the c axis. A number of contacts 
were measured by touching of the film edge after break- 
ing AI2O3 substrate. By this means, the current flows 
preferably along the ab plane. The differential resis- 
tance dV/dl vs V was recorded using a standard lock- 
in technique. The normal resistance Rn (at V ^> A) 
of investigated contacts ranged mainly between 10 and 
1000 n at 4.2 K. 

The important characteristic of PC is their size or 
diameter d, which can be determined from the simple 
formula derived by Wexler P2j for contact resistance: 



Rpc(T) 



16pZ £(T) 



(1) 



where two terms represent ballistic Sharvin and dif- 
fusive Maxwell resistance, correspondingly. Here pi — 
pp/ne 2 , where pp is the Fermi momentum and n is 
the density of charge carriers. The latter for MgB2 
is estimated at n ~ 6.7 x 10 22 [^3j, which results in 
pi ~ 2xl0" 12 O-cm 2 using ^ 5xl0 7 cm/s |§. Hence, 
the upper limit for elastic mean free path I = pl/po for 
our film is about 3 nm. In this case, according to Eq. (1), 
the condition d < I is fulfilled for PC with R > 40 O or 
for lower resistance supposing multiple contacts in par- 
allel. 

We have utilized generally used Blonder-Tinkham- 
Klapwijk equations |24j describing I —V characteristic 
of ballistic N-c-S metallic junctions (here N is normal 
metal, c is constriction and S is superconductor) by ac- 
counting for the processes of Andreev reflection. At fi- 
nite barrier strength at the N-S interface characterized 
by parameter Z ^ and T <C T c , the theory gives the 
dV/dl curves with minima at V ~ ±A/e. To get the 
correct A, the fit of the measured curves to the the- 
ory should be done. The additional smearing of dV/dl 
curves due to, e. g., broadening of the quasiparticle 
DOS in the superconductor can be taken into account 
by including parameter T ]2q| . 

In the case of curves with double gap structure we 
calculated, according to the theory Q , the sum of two 
differential conductances dl/dV with the weight w for 
the larger gap and, correspondingly, with (1 — w) for 
the smaller one. After this, we have transformed dl/dV 
into dV/dl to compare with the measured dependences. 
The best fit was achieved, as a rule, by using its own val- 
ues of Z and T for large and small gap. It is acceptable 
if we suppose that we have a number of microconstric- 
tions with various Z in the region of mechanical touch. 



It is worthy to note that, with increasing of weight fac- 
tor, the difference between Z and T values for large and 
small gap becomes smaller or even vanishes for some 
PCs. 

Results and discussion. Approximately one half (44 
of total 91) of analyzed raw dV/dl vs. V curves show 
visible two-gap structure, although, in most cases with 
shallow features corresponding to a larger gap. The 
samples of some dV/dl curves taken at 4.2 K<CT C with 
double-gap structure, along with calculated curves, are 
shown in Fig. 1. In spite of a number of fitting param- 
eters (A, r, Z, w ) for curves with pronounced (or at 
least visible as shown in Fig. 1) double gap features de- 
termined Al and As are robust as to fitting procedure. 

It turns out that histogram of gaps distribution built 
on the basis of fitting of 44 spectra (see Fig. 2a) has two 
well-separated and quite narrow (especially for the small 
gap) maxima. 

The double-gap distribution is in qualitative agree- 
ment with the distribution of gap values over the Fermi 
surface recently calculated in ( see Fig. 2c). The 
main difference is that theoretical distribution for lower 
gap is wider and has a dominant maxima around 1.6 
meV. This discrepancy can be resolved when we draw 
attention that we have measured curves with double gap 
structure for contacts that is predominantly along the 



c-axis. In this case, according to |26|, gap values along 
c-axis spreads between 2 and 3 meV. The c-axis direc- 
tionality of our measurements is, apparently, the main 
reason of a shallow large gap structure in dV/dl, be- 
cause large gap dominates in the "a-b" plane p6fl . 

It should be mentioned that two very different order 
parameters exist only in the clean limit I 3> 27r£. Since 
in our case I has upper limit in 3 nm and the coher- 
ence length £ ^5 nm [^7j , the observation of two gaps is 
in line with our supposition that in the PC area there 
are small grains with a much larger mean free path. In- 
deed, SEM image of MgB 2 films Jl6| shows that the film 
is granular with 100-500 nm large grains. Therefore, in 
the area of mechanical contact there are some amount 
of small metallic bridges, perhaps, with slightly different 
crystallographic orientation being in parallel. 

The single gap A is seen for the dirty limit and 
is average of small and large gaps with some weights. 
If we assume that this weight has some relation to the 
weight w used in the fitting procedure, then, admit- 
tedly, A ~ w Al + (1 — w)A,s=3.4meV by using upper 
limit w ~0.2 (see Fig. 3). This agrees with the position 
of the maximum of single-gap distribution at 3.5 meV 



3 'ln the case of interface scattering Sharvin resistance should 
be multiple by factor (1+Z 2 )Q. 



4 'On dV/dl of "edge" contacts (a total of 11 curves) only single 
gap structure was observed, probably due to the deterioration of 
the film edge after breaking. 
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(see. Fig. 2b). By the way, according to the calculation 
in Ea ] , a large amount of impurity scattering will cause 
the gaps to converge to A ~4.1 meV. 

Therefore the superconducting properties of this 
compound can be strongly influenced by nonmagnetic 
defects and impurities, which seem to have a great im- 
pact also on the scattering of gap value(s) given by dif- 
ferent authors. 

As to w factor it is hardly to see in Fig. 3 its de- 
pendence on Rn or PC size, which one would expect if 
small gap reflects a degraded surface or large gap is a 
result of surface states Q . 

Table II shows double gap values given by different 
authors. A quite good correspondence between our re- 
sults and data of other authors carried out on different 
types of MgB 2 samples is evident. In our case averaged 
over 44 curves, the ratio of the larger gap to the lower 
one 2.85±0.15 is close to the theoretical value 3:1 fl9fl . 

In conclusion. We have analyzed dV/dl point- 
contact spectra of MgB2 with clear single- and double- 
gap structure. The observed distinct maxima in the 
double gap distribution which is consistent with theo- 
retical calculations |^6| ruled out surface or multiphase 
origin of gap structure and testify about intrinsic su- 
perconducting double-gap state in MgB2- The averaged 
gap value ratio turned out to be in accordance with the 
theoretically predicted ratio 1:3 fl9fl . 
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Figure 1. 

Reduced differential resistance Ti^dV/d/ vs. V mea- 
sured at T = 4.2 K for four MgB2-Ag contacts with 
double gap structure (symbols) . Thin lines are theoret- 
ical dependences calculated with parameters given in 
the Table I. The curves (1-3) are vertically offset for 
clarity. Vertical dashed lines show approximately po- 
sition of large A_l and small As gaps. Experimental 
curves are taken nominally in c-directions 



Table 1. Fitting 


; parameters for curves 


presented 


in Fig. 1. 


Parameters 


Curve 1 


Curve 2 


Curve3 


Curve 4 


Rjv, 


47 


35 


20 


34 


Al, mcV 


7.4 


6.25 


7.35 


7.3 


As, meV 


2.6 


2.54 


2.4 


2.6 


w-factor 


0.11 


0.08 


0.07 


0.06 


Zl 


0.7 


0.71 


0.63 


0.21 


Zs 


0.75 


0.55 


0.56 


0.76 


T L , mcV 


0.4 


0.1 


0.55 





T s , meV 


0.5 


0.54 


0.38 


0.3 



Table 2. Gap values in MgB2 measured by point- 
contact (PCS) or tunneling spectroscopy (TS). 



Method 


Sample 


As, meV 


Ai, meV 


Refs. 


PCS 


Film 


2.45 ±0.15 


7.0 ±0.4 


This pap 


PCS 


Film 


2.3 ±0.3 


6.2 ±0.7 


@ 


PCS 


Grain 


2.8 


7 


(1 


TS 


Granular 


3.9 


7.5 


§ 


TS 


50/i crys. 


3.8 


7.8 




TS 


polycrys. 


1.75 


8.2 


i 



20 



10 



(a) 



A 



3 2 4 6 8 10 



10 



(b) 



j . . i i i . i i 



2 4 6 8 10 



". .ill III. 



(c) 



2 4 6 

Gap value (meV) 



10 



Figure 2. 

The superconducting energy-gap distribution in c-axis- 
oriented MgB2 thin film in the case: (a) double gap 
and (b) single gap. Thin lines show Gaussian fit with 
maxima at (a) 2.45 and 7 meV and (b) 3.5 meV. The 
histogram window of 0.25 meV for (a) and 0.5 meV for 
(b) is chosen to get the most close to normal (Gaus- 
sian) distribution, (c) Distribution of gap values over 
the Fermi surface calculated in pdl 




Figure 3. 

Dependence of the weight factor w on the point-contact 
resistance 



